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using its coupling to a single two-level system. Our approach requires a minimal number of mea- 
surements on the two-level system for a given measurement precision. We consider the application 
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dot: we show that decreasing the electron spin dephasing due to nuclei and increasing the fidelity 
of nuclear-spin-based quantum memory could be within the reach of present day experiments. 
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I. INTRODUCTION 

Decoherence of quantum systems induced by interac- 
tions with low-frequency reservoirs is endemic in solid- 
state quantum information processing (QIP) (HQ- A 
frequently encountered scenario is the coupling of a two- 
level system (qubit) to a mesoscopic bath of two-level 
systems such as defects or background spins. The mani- 
festly non-Markovian nature of system-reservoir coupling 
in this scenario presents challenges for the description of 
the long term dynamics as well as for fault tolerant quan- 
tum error correction 0, The primary experimental 
signature of a low-frequency reservoir is an unknown but 
slowly changing effective field that can substantially re- 
duce the ability to predict the system dynamics. A possi- 
ble strategy to mitigate this effect is to carry out a quan- 
tum measurement which allows for an estimation of the 
unknown reservoir field by controlled manipulation and 
measurement of the qubit. A precise estimation of the 
field acting on the large Hilbert space of the reservoir re- 
quires, however, many repetitions of the procedure: this 
constitutes a major limitation since in almost all cases of 
interest projective measurements on the qubit are slow 
[j| and in turn will limit the accuracy of the estimation 
that can be achieved before the reservoir field changes. 

In this work, we propose a method for estimating an 
unknown quantum field associated with a mesoscopic 
spin ensemble. By using an incoherent version of the 
quantum phase estimation algorithm, @, Q we show that 
the number of qubit measurements scale linearly with the 
number of significant digits of the estimation. We only 
assume the availability of single qubit operations such as 
preparation of a known qubit state, rotations in the xy- 
plane, and measurement, of which only rotations need 
to be fast. The estimation procedure that we describe 
would suppress the dephasing of the qubit induced by 
the reservoir; indeed, an interaction with the estimated 
field leads to coherent unitary evolution that could be 
used for quantum control of the qubit. If the measure- 
ment of the reservoir observable is sufficiently fast and 
strong, it may in turn suppress the free evolution of the 



reservoir in a way that is reminiscent of a quantum Zeno 
effect. 

After presenting a detailed description of the mea- 
surement procedure and discussing its performance and 
limitations, we focus on a specific application of the 
procedure for the case of a single quantum dot (QD) 
electron spin interacting with the mesoscopic nuclear 
spin ensemble defined by the QD. It is by now well 
known that the major source of decoherence for the 
electron-spin qubits in QDs |8( is the hyperfine interac- 
tion between the spins of the lattice nuclei and the elec- 
tron HESmilllllllllll. A particular feature of the 
hyperfine-related dephasing is the long correlation time 
(t c ) associated with nuclear spins. This enables tech- 
niques such as spin-echo to greatly suppress the dephas- 
ing 0. In 0| it was suggested to measure the nuclear 
field to reduce electron spin decoherence times; precise 
knowledge of the instantaneous value of the field would 
even allow for controlled unitary operations. For exam- 
ple, knowledge of the field in adjacent QDs yields an effec- 
tive field gradient that could be used in recently proposed 
quantum computing approaches with pairs of electron 
spins ^tJ. Moreover, with sufficient control, the collec- 
tive spin of the nuclei in a QD may be used as a hi ghly co- 
herent qubit-implementation in its own right |l8lll9ll2fj| . 



II. PHASE ESTIMATION 

In the following we consider an indirect measure- 
ment scheme in which the system under investigation is 
brought into interaction with a probe spin (a two-level 
system in our case) in a suitably prepared state. Measur- 
ing the probe spin after a given interaction time t yields 
information about the state of the system. We assume 
the mesoscopic system evolves only slowly compared to 
the procedure, and further that the measurement does 
not directly perturb the system. In essence, we are per- 
forming a series of quantum non-demolition (QND) mea- 
surements on the system with the probe spin. 
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We consider an interaction Hamiltonian of the form 
Hint = fiA z $ & (1) 

which lends itself easily to a measurement of the ob- 
servable A z . The QND requirement is satisfied for 
[-Hint, H bath] — > 0. The applicability of H mt in situations 
of physical interest is discussed in Sec. I VII Given this in- 
teraction, the strategy to measure A z is in close analogy 
to the so-called Ramsey interfcromctry approach, which 
we now briefly review. 

For example, an atomic transition has a fixed, scalar 
value for A z which corresponds to the transition fre- 
quency. By measuring A z as well as possible in a given 
time period, the measurement apparatus can be locked 
to the fixed value, as happens in atomic clocks. The 
probe spin S is prepared in a state |+) = (||) + \ i))/V2. 
It will undergo evolution under H according to Ut — 
exjp(—itA z S z ). After an interaction time t, the probe 
spin's state will be 

cos(fit) |+) +isin(f2t) |-) (2) 

where 57 = A z /2 is the precession frequency for the probe 
spin. A measurement of the spin in the |±) basis yields a 
probability cos 2 (Sli) of being in the |+) state. Accumu- 
lating the results of many such measurements allows one 
to estimate the value for fl (and therefore A z ). In gen- 
eral, the best estimate is limited by interaction time: for 
an expected uncertainty in A z of Ao and an appropriate 
choice of t, M measurements with fixed interaction times 
1/An can estimate A z to no better than ~ Aq/^/M (see 
|2ll | and references therein). 

In our scenario, the situation is slightly different in 
that A z is now a quantum variable. For a state |s) in the 
Hilbert space of the system Ti. which is an eigenstate of 
A z with eigenvalue 20 s , the coupling induces oscillations: 

U t \s) |+) = \s) [cos(rU) \+)+i sin(rU) |-)] . (3) 

Thus, the probability to measure the probe spin in state 
|+) given that the system is in a state |s) is p(+\s) — 
cos 2 (fl s t) at time t, providing information about which 
eigenvalue of A z is realized. Comparing Eq. J5J to Eq. 
indicates that the same techniques used in atomic clocks 
(Ramsey interfcromctry) could be used in this scenario to 
measure f2 s and thus project the bath in some eigenstate 
of A z with an eigenvalue of fl s to within the uncertainty 
of the measurement. 

Beyond the Ramsey approach, there are several ways 
to extract this information, which differ in the choice of 
interaction times tj and the subsequent measurements. 
The general results on quantum metrology of p2^ | show, 
however, that the standard Ramsey scheme with fixed 
interaction time t is already optimal in that the scal- 
ing of the final variance with the inverse of the total 
interaction time cannot be improved without using en- 
tangled probe states. Nevertheless, the Ramsey scheme 
will not be the most suitable in all circumstances. For 
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FIG. 1: (Color online) Illustration of the first two steps of 
the measurement procedure. The original distribution P(A Z ), 
with rms width Ao is shown at the top. After the first mea- 
surement, with result tm = 0, 1, the conditional distribution 
(middle plots) reflects the knowledge of the least significant 
bit. The next measurement result rjn-i further reduces the 
distribution (bottom plots). 



example, we have assumed so far that preparation and 
measurement of the probe spin is fast when compared to 
1/Ao. However, in many situations with single quantum 
systems, this assumption is no longer true, and it then 
becomes desirable to minimize the number of prepara- 
tion/measurement steps in the scheme. 



III. THE MEASUREMENT SCHEME 

We now show that by varying the interaction time and 
the final measurements such that each step yields the 
maximum information about Q s , we can obtain the same 
accuracy as standard Ramsey techniques with a simi- 
lar interaction time, but only a logarithmic number of 
probe spin preparations and measurements. As a trivial 
case, if A z had eigenvalues and 1 only, then measur- 
ing the probe in the ±-basis after an interaction time 
ii = it, we find +(— ) with certainty, if the system are in 
an A z — 0(l)-eigenstate; if they were initially in a super- 
position, measuring the probe projects the system to the 
corresponding eigenspaces. We can extend this simple 
example (in the spirit of the quantum phase estimation 
algorithm [a, [7j and its application to the measurement 
of a classical field |23|) to implement an ^-measurement 
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by successively determining the binary digits of the eigen- 
value. We start with the ideal case, then generalize to a 
more realistic scenario. 



A. Ideal case 

If all the eigenvalues of A z are an integer multiple of 
some known number a and bounded by 2 M a, then this 
procedure yields a perfect Admeasurement in M steps: 
let us write all eigenvalues as 2fi s = a2 M Xw=i s/2~ ; . 
The sum we denote by s and also use the notation s = 
O.S1S2 . . . Sm- Starting now with an interaction time t\ — 
~, we have f2 s ii = sm § m °d it. Hence the state of the 
probe electron is flipped if and only if sm = 1- Therefore 
measuring the probe electron in state +(— ) projects the 
nuclei to the subspace of even(odd) multiples of a (see 
Fig. . We denote the result of the first measurement 
by tm = 0(1) if the outcome was "+(—)". All the higher 
digits have no effect on the measurement result since they 
induce rotations by an integer multiple of 7r which have 
no effect on the probabilities p(±\s). 

To measure the higher digits, we reduce the interac- 
tion time by half in each subsequent step: = 2~Hi 
until we reach tjw = — in the final and shortest step. For 
j > 1 the rotation angle fl s tj (mod it) in the jth step 
does not only depend on the jth binary digit of s but 
also on the previous digits (which have already been mea- 
sured, giving results r M +i-i = s M +i-i,l = 1, . .., j - 1). 
The angle fl s tj (mod tt) is given by SM+i-jf + fj with 
<Pj = £ J2i=i r A/+i— «2 /— J , where we have used the results 
ri already obtained. This over-rotation by the angle tpj 
can be taken into account in the choice of the measure- 
ment basis for the jth step: if the jth measurement is 
performed in a rotated basis |±j) that is determined by 
the previous results 77, namely 

:= cos ipj |+> - ismipj |-) , (4a) 
:= simpj |+) +icos(pj \— ) , (4b) 

then the jth measurement yields "+" (rM+i-j = 0) if 
SM+i-j = and "— " {vm+i-j = 1) otherwise. Thus, af- 
ter M measurements we obtain 77 = s;, VZ = 1, . . . ,M 
and have performed a complete measurement of A z 
(where the number M of probe particles used is the small- 
est integer such that 2 M > A z /a). 

Before proceeding, we note that the proposed scheme is 
nothing but an "incoherent" implementation of the quan- 
tum phase estimation algorithm: As originally proposed, 
this algorithm allows measurement of the eigenvalue of a 
unitary U by preparing M qubits (the control-register) 
in the state |+}® M (i.e., the equal superposition of all 
computational basis states |j) , j = 0, . . . , 2 M — 1) and 
performing controlled-?/ 23 gates between the jth qubit 
and an additional register prepared in an eigenstate \s) of 
U with U \s) = e l2ns \s). The controlled-t/ gates let each 
computational basis state acquire a s-dependent phase: 
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|^). Then the inverse quantum Fourier trans- 



formation (QFT) is performed on the control register, 
which is then measured in the computational basis, yield- 
ing the binary digits of s. Performing the QFT is still a 
forbidding task, but not necessary here: the sequence of 
measurements in the rotated basis |±j) described above 
is in fact an implementation of the combination of QFT 
and measurement into one s tep. This was previously sug- 
gested in different contexts p lIMI^ . 



B. Realistic case 

In general, there is no known a such that all eigenval- 
ues s of A z are integer multiples of a. Nevertheless, as 
discussed below, the above procedure can still produce 
a very accurate measurement of A z if sufficiently many 
digits are measured. Now we evaluate the performance 
of the proposed measurement scheme in the realistic case 
of non-integer eigenvalues. Since here we are interested 
in the fundamental limits of the scheme, we will for now 
assume all operations on the probe qubit (state prepara- 
tion, measurement, and timing) to be exact; the effect of 
these imperfection is considered in Sec. Without loss 
of generality, let A and denote the largest and smallest 
eigenvalues of A z , respectively ^(| and choose a = 2 A 
such that the eigenvalues of A z /a are all € [0,1/2]. These 
are the eigenvalues s we measure in the following. 

The function from which all relevant properties of our 
strategy can be calculated is the conditional probability 
Pm(R\s) to obtain (after measuring M electrons) a result 
R = 0.rir 2 . . . tm given that the system was prepared 
in an eigenstate with eigenvalue s. The probability to 
measure R is given by the product of the probabilities 
to measure tm+i-j in the jth step, which is cos 2 (fl s tj 



rM+i-j% ) = cos 2 (tt[s - R]2 M -i). Hence 



M-l 



p M (R\s) = Yl cos 2 (^[s-i?]2 fc ), 



(5) 



k=0 



see also |2jj . This formula can be simplified by repeatedly 
using 2sinx = sin(x/2) cos(x/2) to give 



p M {R\s) = 



sin(2 M 7r[s - R}) 
2 M sin( 7 r[s- J R]) 



(6) 



Assume the nuclei are initially prepared in a state p with 
prior probability p(s) to find them in the eigenspace be- 
longing to the eigenvalue s. After the measurement, we 
can update this distribution given our measurement re- 
sult. We obtain, according to Bayes' formula: 



p M {s\R) 



p M {R\s)p{s) 



T, s Pm(R\s)p(s)' 
with expectation value denoted by s/j. 



(7) 
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IV. PERFORMANCE OF THE SCHEME 

As the figure of merit for the performance of the mea- 
surement scheme we take the improvement of the average 
uncertainty in A z of the updated distribution 



AA, 



-(M) 



R 



Pm{R)JY(s~- S b) 2 Pm{s\R) (8) 



over the initial uncertainty Ao = AAz . An upper 
bound to AA ZMS ^ MS> is given by the square root of the 
average variance V 



V M :=^PM(^)^(s-Sfl) 2 PM(s|-R), 



(9) 



as easily checked by the Cauchy-Schwarz inequality. We 
now show that V m < 2~ M . We replace s R — > R = 
min{i?, 1— R}; we can use any such replacement to obtain 
an upper bound, as the expectation value x — ^2 p{x)x 
minimizes v{y) — ^2 x p(x)(x — y) 2 . This choice means 
that measurement results R > 1/2 are interpreted as 
1 — R, which is appropriate since the scheme does not 
distinguish the numbers s = S and s' = 1 — 8 and due to 
the choice of a only s £ [0, 1/2] occur. Thus 

V M < ^p(s)^( S -fc(4)= : E^( s )' 



The terms Vm(s) can be shown |47j to be < by2~ M 
with by = (l + 2~ M )/2. This means that perform- 
ing M measurements yields a state with ^-uncertainty 
AA2. cs t < a2~ M / 2 . For example, we need about 13 in- 
teractions with the probe spin to reach the 1%-level in 
AA Zi e St /a and about 7 more for every additional factor 
of 10. 

The overall procedure requires a total time Tm = 
r m ) = 2ti (1 — 2~ M ) + Mr m , which is an 



interaction time (determined mainly by the time t\ = 
f 2 M = 7r(2/A ) -1 2 M needed for the least significant 
digit probed) and the time to make M measurements 
(r m is the time to make a single measurement). We ob- 
tain for the average uncertainty an upper bound in terms 
of the interaction time T; n t = Tm — Mr m needed: 



AA 2 



< 



AnTj 



O-tint 



(10) 



Immediately the similarity with standard atomic clock 
approaches is apparent, as the uncertainty decreases with 
the square root of the interaction time. However, while 
for an atomic clock scheme, in which the interaction time 
per measurement is kept fixed to ~ 7r/(/Ao), the total 
time to reach the precision of Eq. (|10|) is T M = 2t\ + 
T m 2 M . For our method the measurement time is reduced 
dramatically by a time T T M - T M = r m (2 M - M). In 
this manner our approach requires a polynomial, rather 
than exponential, number of measurements for a given 



binaiy digits probed 

FIG. 2: (Color online) Estimation in the presence of prepara- 
tion and measurement error p: the logarithm of the improve- 
ment Aj4 z ,est/Ao is plotted versus the number M of binary 
digits measured. The solid lines represent different error rates 
(p = 0: black circles; p — 10~ 4 blue crosses; p = 10~ 2 : red 
triangles; and p = 10 _1 : green diamonds). 
The broken red curves (triangles) demonstrate the benefit of 
simple error correction (for p = 10 -2 ): strategy I (3 repeti- 
tions per digit, use majority result; either for all digits or only 
for leading M/2 digits): (dash-dotted, dotted - almost undis- 
tinguishable) ; strategy II (increase number of repetitions for 
more significant digits to 7 repetitions for leading M/8 digits, 
5 for next M/8, and 3 for next Af/4): dashed. We see that 
the latter provides a better accuracy than the uncorrected 
p — 1CP 4 curve. 



accuracy, though the overall interaction time is the same 
for both techniques. 

It may be remarked that even the scaling in interaction 
time differs significantly if other figures of merit are con- 
sidered. For example, our scheme provides a square-root 
speed-up in T; nt over the standard Ramsey scheme if the 
aim is to maximize the information gain or to minimize 
the confidence interval 1271. 



V. ERRORS AND FLUCTUATIONS IN A z 

Up until now we have considered an idealized situation 
in which the value of A z does not change over the course 
of the measurement and in which preparation and mea- 
surement of the probe system work with unit fidelity. Let 
us now investigate the robustness of our scheme in the 
presence of these errors. 



A. Preparation and Measurement Errors 

By relying upon a small number of measurements, the 
scheme we described becomes more susceptible to prepa- 
ration and measurement errors. An error in the determi- 
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nation of the fcth digit leads to an increase of the error 
probability in the subsequent digits. This error amplifi- 
cation leads to a scaling of the final error of y/p, where 
p is the probability of incurring a preparation or mea- 
surement error in a single step. Wc confirm this with 
Monte Carlo simulations of the measurement procedure 
(Fig. EH), leading to an asymptotic bound: 

AA z , cst < (2/A )VP- ( n ) 

Standard error correction (EC) techniques can be used to 
overcome this problem. E.g., by performing three mea- 
surements for each digit and using majority vote, the 
effective probability of error can be reduced to ~ 3p 2 - 
at the expense of tripling the interaction time and num- 
ber of measurements. While this may look like a big 
overhead, it should be noted that the scheme can be sig- 
nificantly improved: the least significant digits do not 
require any EC. For them, the scheme gives noisy results 
even for error rate p = due to the undetermined digits 
of A z , this does not affect the most significant M/2 dig- 
its. This indicates that it may be enough to apply EC 
for the leading digits. 

As can be seen from Fig. Eb, this simple EC strat- 
egy provides a significant improvement in the asymp- 
totic AA Z . This is hardly changed, when EC is applied 
only to the leading half of the digits. Thus only twice as 
many measurements (and an additional interaction time 
~ 2 M / 2+2 which is <C Tm) is needed for an order-of- 
magnitude improvement in AA Z . By repeating the mea- 
surement of more important digits even more often, the 
effect of technical errors can be reduced even further, 
as confirmed by Monte Carlo simulations (Fig. Et>) ■ We 



also note that further improvements (beyond M digits) 
can be achieved by this technique. In essence, choos- 
ing a digital approach to error correction for our digital 
technique yields substantially better performance than 
adapting the digital technique to an analog approach. 



B. Estimation of bath decorrelation errors 

In practice, internal bath dynamics will lead to fluctu- 
ations in A z , such that A z (t) ^ A z (t') for times t and 
t' that are sufficiently different. Furthermore, apparatus 
errors, as outlined above, lead to errors in our measure- 
ment procedure. We will assume that the variations of 
A z are slow over short time intervals, allowing us to ap- 
proximate the M bit measurement process as a continu- 
ous measurement over the time Tm with some additional 
noise with variance ps AA Zyest (M) 2 <C Aq. Then we will 
find the expected difference in our measurement result 
and the value of A z at a later time. 

Under the above approximations the value of the fcth 
such measurement (where a complete set of M bits takes 
a time Tm and the fcth such measurement ends at time 
tk) is 

m k = -L / A z (t)dt + G k (12) 

J-M Jt k -T M 

where the noise from measurement is incorporated 
in the stochastic noise variable Gk with (GkGk') ~ 
8kk' AA z ^ st {M) 2 . We can estimate A z at a later time, 
and find the variance of this estimate from the actual 
value: 



V M (t>t k ) = ([m k - A z {t)} 2 ) 



= (G 2 ) + A 2 + — - / {{Az (t'),A z (t")} + )dt"dt'-— ({A z (t'),A z (t)} + )dt' 

L1 M Jt k -Tu Jtk-Tu ±M Jtk-TM 

If we assume A z is a Gaussian variable with zero mean, described by a spectral function S(u>) (i.e., (A z (t)A z (t+T)) = 
f^Siuiy^duj), then 

1 MJ-oo ( uj / Z ) J-oo J-MU/Z 



For A z that fluctuates slowly in time and corresponds to a 
non-Mar kovian, low frequency noise, the second moment 
of S(ui) converges. We define: 



1 

7! 



i 



S(u>)u> 2 div 



(13) 



When T M ,t — t k + T M -C t c , we may expand the sine 
terms in the integrals. Taking t = t k + Tm, the expected 



variance to order (Tm/^c) 2 is 



V M (t)~ AA z . cst (M) 2 + Al 



'T M ' 


2 


(1 




. ^ _ 




V3 





(14) 



As an example case, we consider as realistic parame- 
ters t c ps 1ms and Tm — 16 ^s with AA 2 zcst /AQ = 
0.025 2 . These parameter choices are described in detail 
in Sec. I VII We find that our variance 16/is after the mea- 



6 



surement is approximately 0.035 2 Ag with equal contri- 
butions from the measurement noise and from the bath 
decorrelation. Substantially faster decorrelation would 
dominate the noise in the estimate, and render our tech- 
nique unusable. 

In the limit of slow decorrelation, this approach would 
allow one to use the (random) field A z to perform a con- 
trolled unitary of the form exp(— irrikS z T) at a time t, 
with a fidelity 

F = exp(-{(/ A z {t')dt' -m k f)/A) (15) 

Jt-T 

For example, a 7r rotation around the probe spins' z axis 
would have a fidelity «1- V M (t)n 2 / A 2 ,, or 0.998 for the 
above parameters. 

We remark that this approach for estimation in the 
presence of bath fluctuations is not optimal (Kalman fil- 
tering ,28J would be more appropriate for making an es- 
timation of A z using the measurement results). Further- 
more, it does not account for the non-linear aspects of 
our measurement procedure, nor does it incorporate any 
effect of the measurement on the evolution of the bath 
(e.g., quantum Zeno effect). More detailed investigations 
of these aspects of the process should be considered in an 
optimal control setting. Nonetheless, our simple analysis 
above indicates that slow decorrelation of the bath will 
lead to modest additional error in the estimate of A z . 



VI. EXAMPLE: ESTIMATING COLLECTIVE 
NUCLEAR SPIN IN A QUANTUM DOT 

Now we apply these general results to the problem of 
estimating the collective spin of the lattice nuclei in a 
QD. 

The interaction of a single electron spin in a QD with 
the spins of the lattice nuclei Ij is described by the Fermi 
contact term |nj 

S-£>4 (16) 

3 

where the sum in Eq. (|16fl runs over all the JV lattice nu- 
clei. The aj are constants describing the coupling of the 
jth nuclear spin with the electron. They are proportional 
to the modulus squared of the electron wave function at 
the location of the jth nucleus and are normalized such 
that J2j OLjlV' = A, which denotes the hyperfine cou- 
pling strength. 

Due to the small size of the nuclear Zeeman energies, 
the nuclei are typically in a highly mixed state even at 
dilution refrigerator temperatures. This implies that the 
electron experiences an effective magnetic field (Over- 
hauser field, J? n uc) with large variance, reducing the fi- 
delity of quantum memory and quantum gates. This re- 
duction arises both from the inhomogencous nature of the 
field (i?nuc varies from dot to dot) [2^| and the variation 



of B nuc over time due to nuclear-spin dynamics (even a 
single electron experiences different field strengths over 
time, implying loss of fidelity due to time-ensemble aver- 
aging). 

In a large external magnetic field in the z-direction the 
spin flips described by the x and y terms are suppressed 
and - in the interaction picture and the rotating wave 
approximation the relevant Hamiltonian is of the type 
given in Eq. (JIJ, where A z is now the collective nuclear 
spin operator 

JV 

A z = £ a ^ j \ (17) 
i=i 

which gives the projection of the Overhauser field along 
the external field axis by B nuc z = HA z /g* . Before contin- 
uing, let us remark here, that one can expect to obtain an 
effective coupling of the type Eq. QJ in a similar fashion 
as a good approximation to a general spin-environment 
coupling S ■ A, whenever the computational basis states 
of the qubit are non-degenerate (as guaranteed in the sys- 
tem studied here by the external field) and the coupling 
to the environment is sufficiently weak such that bit-flip 
errors are detuned. 

To realize the single-spin operations needed for our 
protocol - preparation, rotation, and read-out - many 
approaches have been suggested as part of a quantum 
computing implementation with electron spin qubits in 
QDs using either electrical or optical control (see, e.g., 
30] for a recent review). 

The experimental progress towards coherent single spin 
manipulation has been remarkable in recent years. In 
particular, the kind of operations needed for our proto- 
col have already been implemented in different settings: 
For self-assembled dots, state preparation with F > 0.99 
has been realized [3l|], while for electrically defined dots, 
single-spin measurement with a fidelity of F > 0.72 was 
reported In the double-dot setting all three 

operations have recently been demonstrated, and we es- 
timate the combined fidelity to be F > 0.7. 

As can be seen from Fig. El at the level of 1% accuracy 
of state preparation, rotation and read-out, the proposed 
nuclear spin measurement should be realizable. As dis- 
cussed in many specific proposals [3(j these error rates 
appear attainable in both the transport and the optical 
setting. Apart from single qubit operations, our proposal 
also requires precise control of the interaction time. Fast 
arbitrary wave form generators used in the double-dot 
experiments, have time resolutions better than 30 ps|4S| 
and minimum step sizes of 200 ps, which translates into 
errors of a few percent in estimating A z with initial un- 
certainties of order 1 ns _1 . Uncertainties of this order are 
expected for large QDs (JV ~ 10 6 ) even if they are un- 
polarizcd and for smaller ones at correspondingly higher 
polarization (see below). 

For GaAs and InAs QDs in the single electron regime, 
A ~ 50 - 200ns" 1 and JV ~ 10 4 — 10 6 . The uncertainty 
A 2 , = (Al) - (A z ) 2 = (T 2 *)~ 2 determines the inhomoge- 
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neous dephasing time T 2 * |12( • Especially at low polariza- 
tion P, this uncertainty is large Aq sw A 2 (l — P 2 )/N, and 
without correction iJ zz leads to fast inhomogeneous de- 
phasing of electron spin qubits: T 2 * ^ 10 ns nas been 
observed Q |U 13 • However, as A z is slowly vary- 
ing 0, 0| , it may be estimated, thereby reducing the 
uncertainty in its value and the corresponding dephas- 
ing. This is expected to be particularly effective, when 
combining estimation with recent progress in polarizing 
the nuclear spin ensemble |33l l35l l3t| . 

In a QD system such as with r m ~ 1/is and for 
1/Aq — 10ns we can estimate 8 digits (M = 8) (improv- 
ing A^4 ZjOSt by a factor of at least 16) in a total time 
Tm = 16/xs. In contrast, a standard atomic clock mea- 
surement scheme would require a time ~ 280/iS. 

We now consider limits to the estimation process, fo- 
cusing on expected variations of A z due to nuclear spin 
exchange and preparation and measurement errors. Nu- 
clear spin exchange, in which two nuclei switch spin 
states, may occur directly by dipole-dipole interactions 
or indirectly via virtual electron spin flips. Such flips lead 
to variations of A z as spins i and j may have a j ^ ctj . 

The dipole-dipole process, with a 1/r 3 scaling, may 
be approximated by a diffusive process at leng th scales 
substantially longer than the lattice spacing |l5ll37j . The 
length scale for a spin at site i to a site j such that a.j ~ 
on is not satisfied is on the order of the QD radius (5-50 
nm); for diffusion constants appropriate for GaAs |3Sj . 
the time scale for a change of A z comparable to AA Z by 
this process is ~ 0.01 — 10s. 

However, nuclear spin exchange mediated by virtual 
electron spin flips may be faster. This process is the 
first correction to the rotating wave approximation, and 
is due to the (heretofore neglected) terms in the contact 
interaction, H$ = ^(A + S- + A-S+), which are sup- 
pressed to first order by the electron Larmor precession 
frequency e 7 ,. These have been considered in detail else- 
where [IS 0111 EE! El- Using P erturbation theory 
to fourth order, the estimated decorrelation time for A z 
is t^ 1 — A 2 /(e z N 3 / 2 ), giving values 0.1 — 100ms -1 for our 
parameter range Taking t c = 1 ms, we may estimate 
the optimal number of digits to measure. Using Eq. (|14f> . 
the best measurement time is given by Tm ~ t c /2 M / 2 and 
for the values used above, M m 10 — 11 is optimal. We 
note as a direct corollary that our measurement scheme 
provides a sensitive probe of the nuclear spin dynamics 
on nanometer length scales. 

We now consider implications of these results for im- 
proving the performance of nuclear spin ensembles, both 
as quantum memory [l9l | and as a qubit |20j . The domi- 
nant error mechanism is the same as for other spin-qubit 
schemes in QDs: uncertainty in A z . The proposed mea- 
surement scheme alleviates this problem. However, the 
nuclear spin ensembles operate in a subspace of collec- 
tive states |0) and |1), where the first is a "dark state", 
characterized by A^ |0) = (and the second is cx A + |0), 
where A± — ajl±). Thus |0) is an A^ eigenstate 
and cannot be an A z eigenstate when ak ^ const (ex- 



cept for full polarization). Therefore, the measurement 
[which essentially projects to certain A z -"eigenspaces" 
(A z £ [a — Aa, a + Aa])] moves the system out of the 
computational space, leading to leakage errors. The in- 
commensurate requirements of measuring A z and using 
an A- eigenstate place a additional restriction on the 
precision of the measurement. The optimal number of 
digits can be estimated in perturbation theory, u sing an 
interaction time T mt ~ 2ti and numerical results |l9j on 
the polarization dependence of A . We find that for high 
polarization P > 90% a relative error of Aa/a < 1% is 
achievable. 



VII. CONCLUSIONS 

We have shown that a measurement approach based 
on quantum phase estimation can accurately measure 
a slowly varying mesoscopic environment coupled to a 
qubit via a pure dephasing Hamiltonian. By letting a 
qubit interact for a sequence of well controlled times and 
measuring its state after the interaction, the value of the 
dephasing variable can be determined, thus reducing sig- 
nificantly the dephasing rate. 

The procedure requires fast single qubit rotations, 
but can tolerate realistically slow qubit measurements, 
since the phase estimation approach minimizes the num- 
ber of measurements. Limitations due to measurement 
and preparation errors may be overcome by combining 
our approach with standard error correction techniques. 
Fluctuations in the environment can also be tolerated, 
and our measurement still provides the basis for a good 
estimate, if the decorrelation time of the environment is 
not too short. 

In view of the implementation of our scheme, we have 
considered the hyperfine coupling of an electron spin in 
a quantum dot to the nuclear spin ensemble Our cal- 
culations show that the Overhauser field in a quantum 
dot can be accurately measured in times shorter that the 
nuclear decorrelation time by shuttling suitably prepared 
electrons through the dot. Given recent advances in elec- 
tron measurement and control 0, llGj this protocol could 
be used to alleviate the effect of hyperfine decoherence of 
electron spin qubits and allow for detailed study of the 
nuclear spin dynamics in quantum dots. Our approach 
complements other approaches to measuring the Over- 
hauser field in a quantum dot that have recently been 
explored [illi^. 

While we discussed a single electron in a single quan- 
tum dot, the method can also be applied, with modifica- 
tion to preparation and measurement procedures |49l , to 
the case of two electrons in a double dot Q, l34l |45|. 

As we have seen, the Hamiltonian Eq. can serve as 
a good approximation to more general qubit-environment 
coupling in the case of weak coupling and a non- 
degenerate qubit. Therefore, we expect that this tech- 
nique may find application in other systems with long 
measurement times and slowly varying mesoscopic envi- 
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ronmcnts. 
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